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Abstract 

In this paper we first derive a Coulomb Hamiltonian for electron- 
electron interaction in quantum dots in the Heisenberg picture. Then we 
use this Hamiltonian to enhance a Bloch model, which happens to be 
nonlinear in the density matrix. The coupling with Maxwell equations 
when interaction with an electromagnetic field is also considered from the 
Cauchy problem point of view. The study is completed by numerical 
results and a discussion about the advisability of neglecting intra-band 
coherences, as is done in part of the literature. 

Keywords: Maxwell-Bloch model, quantum dot. Coulomb interaction, Cauchy 
problem, Liouville model, positiveness properties. 

1 Introduction 

Bloch model is a very common model to describe the time evolution of a system 
of electrons in different contexts such as gases of electrons, glasses or crys- 
tals. The very classical case of gases and glasses involves isotropic media. The 
electrons are supposed to be localized and non interacting. Their behavior is 
averaged at the mesoscopic scale. This leads to relatively simple models where 
matter energy levels are quantized and labelled by integers. The case of crys- 
tals |BBFB-|-04l also involves levels integer indexed levels, but symmetries and 
directions in matter have to be taken into account. 

Bloch model has also been extended for the description of quantum wells 
|HK96[ IKR92j . and quantum dots |GH02[ IBFIO] , In these models matter is 
described by the state of two species of particles (electrons and holes, or equiv- 
alently valence and conduction electrons). In quantum wells, energy levels are 
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indexed by vectors. In contrast the confinement of electrons in quantum dots 
leads to integer indexed levels like in gases, which often leads to consider quan- 
tum dots as pseudo-atoms, but this is a very raw vision. In particular, among 
other differences, electrons at the same mesoscopic location do interact directly 
via Coulomb interaction. 

[BFlOj is a preliminary paper that derives raw Bloch equations for two 
species of electrons (conduction and valence) only taking into account the free 
electron Hamiltonian and the interaction with a laser electric field. The aim 
of the present paper is to include properly Coulomb interaction in this model. 
Beyond the sole derivation of the model, we want in particular to study its math- 
ematical properties. In the continuation of |Bid01[ |BBR01| . we want to show 
that a certain number of properties are preserved through the time evolution, 
such as Hermicity and positiveness of the density matrix. 

1.1 The raw Bloch model 

Let us first recall the main results obtained in [BFlOj and fix the notations. 

1.1.1 Commutators and Heisenberg equation 

Let A and B be two operators, we define their commutator by [A, B] = AB—BA 
and their skew-commutator by {A, B} = AB ~\- BA. For an operator A, we 
will define the associated observable {A) = Tt{SqA) by averaging with respect 
to the initial state density 5*0 of the system. If the system is described by an 
Hamiltonian iJ, the time-evolution for this observable is given by the Heisenberg 
equation 

mdt{A)^{[A,H]). (1) 

When the observable is the density matrix, the Heisenberg equation of motion 
is called the Bloch equation. 

1.1.2 Operators for quantum dots 

A quantum dot is defined as a collection of conduction and valence electrons. 
There is of course no conduction in quantum dots since the electrons are con- 
fined in every direction, but this terminology is useful to distinguish between 
the valence electrons — which are in fact the absence of holes in the valence 
band, see Section [2?2| — and the free, but confined, electrons. For each species, 
energy levels are quantized and indexed by a set of integers, /"^ and P , for con- 
duction and valence electrons respectively. For i € P , we define the creation 
and annihilation operators cj and Electrons are fermions and should respect 
the Pauli exception principle. The corresponding skew-commutation rules are 

{c,;,Cj} = {4,c]} = 0, {Cr,c]} = 6ij. 

This implies in particular that CiCi — clcl — 0, which means clearly that it is 
impossible to create twice or annihilate twice the same electron. This is the 
Pauli exclusion rule. 
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Likewise, for valence electrons, we define the creation and annihilation op- 
erators v] and Vi, for which 

{v^ ,vj} = {vl,v]} = 0, {v, ,v]} = (5, J . 

Of course any conduction operator commutes with any valence operator. 

1.1.3 Observables for quantum dots 

The observable we are interested in is the density matrix. It includes a conduc- 
tion density matrix, which elements are the p^^ = {cjCi). This matrix is clearly 

Hermitian. Its diagonal terms — (c|^Ci) are also called populations and give 
the probability to find an electron in state i. The off-diagonal terms, p^^ , i j 
are called (intra-band) coherences. Of course we also define a valence density 
matrix, which elements are the plj = (ujui). Besides we are interested in inter- 
hand coherences defined by pJJ = ("yjci). The entries of these matrices are the 
variables of the Bloch equation. They are cast in a single density matrix 

P \ ) 
where p'^'^ = p'^^* , which ensures that p is Hermitian. 

1.1.4 Free electron Hamiltonians and interaction with a laser 

The raw Bloch equation for quantum dot is derived in [BFlOj . It takes into 
account two types of Hamiltonians in the Heisenberg equation for the density 
matrix, namely free electron Hamiltonians and interaction Hamiltonians with a 
laser field. The free electron Hamiltonians read 

kel" kGl^ 

for conduction and valence electrons respectively. The coefficients and are 
the free electron energies associated to each electron level (eigen-state) . The 
interaction with a laser characterized by its time-dependent electric field E(f) 
is described by the Hamiltonians 

H"^' = I E (EW-M^,4Q+E*(i).M^,tc|c,), 

H"-" ^ I E m)-Ml,vlvi+B*{t)-Mi:vlv,), 
{k,i)eii-r 

^Lcv ^ ^ (E(t).M-4^;,-|-E*(i)-M-*^;t^,), 

(fc,/)e/'=x/^ 
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where the dipolar moment matrices may be expressed in terms of the wave 
functions associated to each leveh 

Mil = ldr^rir)rrk{r), 

Injecting these Hamiltonians in the Heisenberg equation, the raw Bloch equa- 
tions can be cast in Liouville form 

mdtp^[v^{t),p]. (2) 

In equation Fo(t) = V^ + V^{t) is a sum of a constant term stemming 
from the free energies collected in diagonal matrices Eq = diag({ej}ig7c) and 
-Eg = diag({ej'}ig7v), and a time dependent term due to the interaction with 
the electric field: 

V^-( M and F^fi) - ( ■ ■ ^ 

1.1.5 Mathematical properties of the Liouville equation 

Equation ^ clearly preserves the Hermitian structure of p. Its exact solution 
is ^ ^ 

p{t) - exp (^-^ Fo(t) dr^ p(0) exp Q ^ Fo(t) drj . (3) 

This expression allows to prove a certain number of properties (see [BFIO] ). in 
particular: 

• Equation (|2| is globally well-posed, 

• for all time p{t) is a positive matrix, 

• its trace is conserved through the time evolution. 

1.2 Issues related to the Coulomb Hamiltonian 

The outline of this paper is as follows. In Section [2] we derive the Coulomb 
Hamiltonian in terms of the conduction and valence operators. The associated 
Heisenberg equation is derived in Section [3] but it ends up with an open system 
of equations. The system is closed in Section ?? using Wick theorem, and the 
final Bloch equation has a Liouville form, but is nonlinear. This nonlinearity 
does not allow to use previous literature directly and we prove anew the Her- 
micity, positiveness and boundedness results in Section [4] In Section [s] we show 
the impact of the Coulomb contribution in numerical results and also compare 
our model with a vanishing inter-band coherence model defined in |GH02j . 
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2 Second quantification Coulomb Hamiltonian 



2.1 Electron— hole model 

Coulomb interaction is easier to introduce first in the electron-hole picture and 
using field operators. We denote by 'ipl{r) and tpl^{r) the creation field-operators 
of respectively an electron and a hole at the space location r, and ijje{r) and 
i^hir) the corresponding annihilation field-operators. Then considering that 
there are N relevant electrons in the quantum dot, we can simply write the 
Coulomb Hamiltonians as 

N 

2 



= \T. 11 d'^^dr, p,{n)^l{ri)V''{n,r,)Urj)Mri), (4a) 

N 

= 2 E // dndr,- ^i{n)i>i{ri)V\ru r,)^h(r,)4(n), (4b) 



where V^, and V"-^ are the electron-electron, hole-hole and electron-hole 
Coulomb potentials. The total Coulomb Hamiltonian is H'^ = H^^^ + H^^^^ + 
/fe-h want to derive Bloch-type equations for the Coulomb interaction. 
Bloch equations have the advantage not to depend explicitly on the exact form 
of the field-operators. To this aim we write 

Mr) = E V^Kr-) = Yl CWcL (5a) 

Mr) = E Vil(r-) = Yl ^l*i'')dl (5b) 

In this decomposition (rcsp. da) is the annihilation operator for electrons 
(resp. holes) and cj, (resp. rfj,) is the corresponding creation operator, where 
a take its values in the set of indices P (resp. I^) for electrons (resp. holes) 
states. They are weighted by the electron and hole wave functions Va(^) 
'ip\{r). In the sequel, to avoid unnecessary written complexity, we will often 
omit to specify which set the indices belong to. Pauli exclusion principle is 
ensured by skew-commutation properties. Of course electron operators commute 
with hole operators. The properties are expressed using skew-commutators 
{A, B] = AB + BA and the Kronecker symbol 5. 

Property 1. 

{Ca,c\} = 5ap, {c„,C^} = {4,C^} = 0, 

{doc.di} = 5ai3, {da,dij} = {d\,di} = Q. 
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Inserting decompositions Q in Hamiltonians Q we obtain 



a'i,a'2 



^h-h^ ^ i?:i,.,.i„i4,4.rf.^rfai, (6b) 
a'i,a'2 

H''-'^= Y R'-^ , ,cldld^,c^,, (6c) 



where 



^a.a.aWi = C (^C ('^O^'^ ('^^ '^OCi ('^OCi (0 , (7a) 

^ ?// <W<(^')^'(^'^')<(-')e,(^), (7b) 

^^Ta.aW;. ^ ^'// C(0<(^')^'=~'(r-,r')<(/)Ci(0- (7C) 

The symmetries in the integrands induce the following properties. 

Property 2. Since V{r, r') is an even function of r — r', variables r and r' play 
the same role and 

pc pe ph ph 

Since V{r, r') is a real valued function 

ne _ pc* ph _ ph* 

QiQ2CKiCt2 Q'^Q2'^ia2' aia2CK^a2 a^a2aiQ2' 

^c— h j^c—h^ 

aia2C^iCt2 a'^a2aia2' 

The Pauli exclusion principle (skew-commutation, Property [T]) also induces 
that some terms in i?®"'' and H^~^ are necessarily zero. 

Property 3. If ai — 0L2 or a!^ — a^, 

c 'C I = {) d) d ' d ' = () 

"^01^02^02 "1 ' "'ai"'Q2 "2 "1 

The Coulomb interaction is intended to be included in a global Bloch-type 
model. It has been shown in |BF10| that the electron-hole formulation is not 
proper to write a Bloch model and that a conduction-valence electron model 
should be used instead. 



6 



2.2 Conduction— valence electron model 



The electrons described in the previous section are viewed as conduction elec- 
trons. The presence of a hole in the valence band is the absence of the cor- 
responding valence electron. Creating a hole is annihilating a valence electron 
and vice-versa. We can therefore recover creation and annihilation operators 



for valence electrons already used in Section 1.1.2 as = do, and Va — d\^. 

The order of the operators (which has to be read from the right to the left) 
in the Coulomb Hamiltonians Q has a meaning: in order that two particles 
interact via Coulomb interact they have to pre-exist at locations and rj . Then 
they are annihilated while interacting and recreated at the same locations. The 
transformed Coulomb Hamiltonians ([6]) obey the same rule. 

Definition 1. A product of operators will be said to be in the normal order ^ 
if the annihilation operators are on the right and the creation operators on the 
left. 

If we replace hole operators by valence electron operators in the Hamiltonians 
([6|, we do not obtain normal ordered forms. 

ai ,Q2 : 



Notation 1. In computations we underline the part of the expression which is 
transformed at the next step, e.g. 

For the reader's convenience, many details of the computations using the 
(skew-)commutation rules and Notation [l] are postponed to an appendix. 

To express i?'^^'' we need to compute a normal ordered form of Va-^Va^v^^, w^, 
(see details in Appendix [A|) 

Thanks to Property 2 F^l ^ = ^ , therefore — (5q,q' w^/ Wq, and 
— <5q2Q2^1' ^Qi lead to the same contribution. The same argument can be applied 
to ^a^ai^I^"""! and (5„j„^w^,^Wa2- Hence 
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In the definition of ^ we have dropped the 5a^a'^5a^a'^ and —Sa^a'^^aia'^ 
terms which would lead to zero contributions in the Heisenberg equation. 
In the same way cl^^Va^v^^i^Ca'^ = 5a^a'^c]^^Ca'^ - c^aiV^a'Va2Ca{, hence 

a,Q',/3 ai,a2, 

Setting 

we can define the Coulomb Hamiltonian in the conduction-valence electron 
picture 

OL\,OL2 , 

a'l,a'2 



^^-^= E KZ2o.',.A<^^'.^»i - E ^r,;',4co'. (8c) 



3 Formulation of the Heisenberg equation 

We now write Heisenberg equation ([T]) where A are operators CjCi, VjVi or wjc^ 
and H are the Hamiltonians defined by Equation 

3.1 Wick theorem 

For the free Hamiltonians and the laser interactions the computation of the com- 
mutators led to expressions in terms of the two-operator densities, and therefore 
to a closed set of equations jBFlO] , This will not the case any more here, since 
the commutators stemming from Coulomb Hamiltonians will give rise to four- 
operator densities. To go further we should a priori have evolution equations 
for these observables via the Heisenberg equation, computing commutators with 
the already defined Coulomb Hamiltonians. This would lead inevitably to six- 
operator densities, and so on. To avoid this endless procedure, we have to close 
the system at some point. This is the goal of Wick theorem |Wic50| which 
amounts in our case to write the four-operator densities as sums of products of 
two-operator densities following e.g. the rule 

(WT) 

where the symbol = means " is approximated through Wick theorem by" . 
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3.2 Commutators involving c^Cj 

In order to derive the Heisenberg equation we have to compute [c|ci, H'^~'^] and 
[CjCi, H'^^^] (since [cjCi, H'^"'^] is clearly zero). 

3.2.1 Commutator with 

According to Equation (|8^) 

Remark 1. We already know many situations where [c]ci, cjj^cJj^Cc^Cc'J is nec- 
essarily zero: 

• if none of the indices ai, a2, a'l, is equal either to i or j, 

• if ai = OL2 or a']^ = a2 (see Property [s]), 

We compute separately each commutator (see details in Appendix [b|) 
Using Property [2] we see that 

"i."2 

= 2 'S^ , ,rW / r / - 2 'S^ 7?^= ,r''' r''' 



Ol,Q2, 



We now apply Wick theorem which leads to 

(WT) 



[c^,c.,H^-^]) ^-^^^ 2 5: ((4^";) (^^^'^"2) - {4'^^'2){cM, 

a, 
a[,a'2 

-2 ^aia2ja' ((cliC*) (4.2Ca') - (4iCa') 

ai,a2, 
a 

a.a' 



A:, 
a,a' 
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Defining matrix A'^(p) as 

(Q,a')e(-f")^ 

we can cast the result as 

3.2.2 Commutator with H"-'' 

According to Equation 

Is 



To evaluate the commutators we compute 

Applying this to (a, a') — {ai,a[) in the first sum we obtain 



a' .,13 a,P 



E ((4-"'.) (-^-i) - (4-"i> (-i^-O) 

a. 



a', 13 a, 13 



E ('^i^.') {c]ck) - E (4^"') 

a, a a, a 

- E ^fcojo' (4«"') (4c.) + E Kl-c' {civa') Uc 



k, k, 
a, a Q:,a 



-E^^(4'^^)+E^^^]; 

k,l3 k 
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Using the matrices C^{p), 7^^ where 



and noticing that 



Qkip) = E R^aL'Pl'a, (10) 

1-riP) = - E (11) 

= - E ^^/sfc"/?' (12) 



Epc— V vc „,^,c— v^-.N* 

^akja'Pa'a — Ikj \P) 



we can rewrite this 

(WT) 

k k 

which will prove useful for the matrix formulation in Section [3.5[ 

3.3 Commutators involving vjvi 
3.3.1 Commutator with i/^^ 

According to Equation ([sjs) , we have to two evaluate two types of commutators 
to evaluate [vjvi, H^^^]. The first commutator is clearly computed in the same 
way as [cjc;, H'^"'^] replacing conduction electron operators by valence ones. For 
the second part, we have to compute [vjvi,vl^Va']. Since 

we have 

a,a' ,f3 

a', 13 a,/3 
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Hence 

a, Qi,a2, 
a[,a'2 a' 

a',/3 a,l3 

where A'^(p) and are defined as 

^Ikip) = 2 ^ (Riaka' - Rlaa'k) Pa'ai (13) 
(Q,Q')e(-f")^ 

l^ik — 2 ^ {RJiikp ^ R'}sif3k)- (14) 



3.3.2 Commutator with H"-'' 

According to Equation (jsj:), the commutator [vjvi, H'^^^] only involves the first 
term of H'^^^ , which we can write as 

ai,Q2, 

Setting (a, a') = {a2,a2) in the previous commutator computation, we deduce 
that 

Q, Qi,a2, 

Q, 

a'i,Q2 

- ^a;a2a'j(«c„-)(t;i^t;,)-(4^«,)(«i^c„.)) 

Q:i,a2, 
a' 

= E (4ca') - E ^oTTa' (i^l 

A;, A;, 
a.a.' a, a' 

- E Kii', Uv^) + E ^fcoo'. (^"^"') 



fc, fc, 
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Anew, we can rewrite this 

(WT) 



k k 

using 7'^~^(p) and the new matrix C^{p) defined by 

Qkifi)^ E KZ'kP'^c'^- (15) 

(Q,Q')e(/=)2 

3.4 Commutators involving f jcj 
3.4.1 Commutator with H^^^ 

According to Equation (|8^) 



"1 ,02 I 



We compute 

CzCq,^ Cq,^ Cq^ ^iai '-'ai '-'i*-'a2 ^c^2 ^c^' 



(^zai *^Q2 ^^^l ^*CK2 '^ai '^ct2 '^ct'i ~^ '^**^ct2 



^iai^]y:^^aL^ai ^ia2 '^ai '-'Cto '-'ck'i ^~ ^ai ^^09 '-'a^ '^ck'i ■ 



Hence, using Property [2] and Wick theorem, 



a, a' 



We recognize 



3.4.2 Commutator with ff^"^ 

In the same way, according to Equation 



Qia2CK-LCK2 '" J ' ^2 



+2 E (^^aC^g - -R^a/3a')[t']>vi^^a']': 
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and 



Hence, using Property [2] 

a[W2 a' 

Since wjwJ^Wc' — — w^wjwa' — Sja'vl^ + w^Uq'w], we have 

a^a' ,P a, (3 

Hence 

ai,a2, a,l3 
a' 

a, (3 

k, 
a, a' 

fc,/3 

We recognize 

fe k 

3.4.3 Commutator with i/'^"^ 
According to Equation (Ist 



L J * ' J / V aia2CK-La2'- J "^i "2 "i-l 



CKl ,a2, 



a,Q' ,/3 
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To evaluate the first sum (see details in Appendix [C|) 
and 



ai ,a2, 



IjRSt, C^C^Cq,' — ^m^Q-' C^^CiCct' — ^ia^a' C^Cq-'C^, Rlld 



Hence 



a\a2a' J Oil oi2 

CKl ,02 , 



a, 

- E ^a;:2o'.((4c.)R^.'>-(4^"'>«^0) 

Qi ,a2, 
a' 



ct',/3 



E ^-fco' (^^-^"') ("]^'^-) - E ^-^'fe (^1 



a, a Q,Q 



E ^'fc^^'. (4c.) + E (4^"') 



k, k, 
Q,a a, a 



k,P 



bjc, ^"^^^ QMPkj+Y. ^':rip)pi,-J2 pikiir^p^Y. pTkCiAp)+Y. ^^k" 



where we recognize 

ffc-vA (WT) 

fc 



pTr 
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3.5 Matrix formulation 

We would like to cast the former results as 



c]v, v]v. 



where 



V^-^p) V^p) 



which implies 

([c]c„i?C]) ^=^^ V^p)p^ + V^-^{p)p^^-p^V%p)-p^^V^-'^{p), 



^ V^-%p)p^'' + V^p)p^ - p-'^V^-^'ip) - P^V^^(p), 
^^^^ V%p)p'^ + V''-^{p)p^ - p'V^-^p) - P^^Vip). 



Identifying the coefficients computed in Sections |3.2[ |3.3[ and |3.4[ we obtain 
that the result can indeed be cast as \V'^ [p) , p] where 

V^p) = A%p) + C{p)+^^-\ 
V'-^p) = I'-^P), 

V^-%p) = j^-^*ip) = v^-^*ip), 

V^p) = A^p) + C{p) + ^\ 

and the various matrices have been defined by Equations ([9| to (15 1. 

The evolution equation for the density matrix including also the free electron 
Hamiltonians, the interaction with a laser field, and Coulomb interaction can 
therefore be cast in a Liouville form 



ihdtp=[v{t,p{t)),pl 

where Vit,p{t)) = Vo{t) + V^ipit)). 
3.6 Energy shifts 

Recall (see Equation ^) the raw Bloch equation 

ihdtp^[Vo{t),p], 

where 



(16) 



Vo{t) = 



+ 5RE(t) • M'^ 
E*(i) • M'^^* 



E{t) ■ M^^ 
+ 3?E(i) • 



The symmetries in the definition of M'^ and imply that their diagonal entries 
are zero. On the contrary for Eq and Eq there are only diagonal entries. In the 



16 



evolution equation of p1^, the term only involves Eq, and the other entries 
of p play a role through M'^ and M^^. 

We can therefore easily analyze the Coulomb contributions in terms of shifts 
on the free electron energies and off-diagonal terms. Hence we compute 

^<h{p) = 2 ^ {Riaia' - Rlaa-i) Pa'a+ ^ ^laia' Pa' a ^ ^ Rifizfi ^ 
{a,a'}<£{I''y (Q,Q')e(-f^)2 ;96/^ 

^<h{p) = 2 ^ {Riaia' - Riaa'i) Pl'a+ ^ ^'aia'iPa' a + ^ i^lu 

We can define the energy shift matrices as 5E'^{p) — diag({^ej(p)}ig/c) and 
5E-{p) - diag({5enP)}»e/v). Hence V{t,p{t)) - E{p{t)) + R{t, p(t)), where 
E{t,p{t))^dmg{E-{t,p{t)),E^t,p{t))) and 

E^{p{t)) = El+ 5E^[p{t)) and E^ {p{t)) = E^ + SE^p{t)). 

3.7 Vanishing intra-band coherence assumption 

The model for quantum dots given in |GH02| has clearly been derived mimicking 
quantum well models as in [KR92j . In quantum wells electrons and holes can 
only interact if they "see each other" long enough, i.e. if they have (and are 
indexed by) the same wave vector. This leads morally to weakly coupled two- 
level systems. Therefore, in |GH02| . the variables are the level populations and 
the inter-band coherences. In comparison with our model there are no intra- 
band coherence p'^j or pj'^ for i ^ j. 

We can therefore wonder what becomes of our model if we set intra-band 
coherences to zero. First we notice that in the evolution equation of e.g. p^^-, 
there is e.g. a contribution of Vfj (p^^ — pl^ ) . Hence, even if intra-band coherences 
are initially zero, there are not always zero through the evolution with Equation 



(16). Setting artificially intra-band coherences to zero therefore destroys the 
Liouville structure of the system. 

In our model, this hypothesis also changes the definition of A^{p), C^(p) and 
A'^(p), which become 

^ikip) = 2 {Rjaka ~ Rjaak) Paa^ 

ael" 

Cikip) = ^ ^iakaPaa^ 
^ikip) = 2^ i^iaka ~ ^iaak) Paa- 

4 Mathematical analysis 

4.1 Estimates on the density matrix 

In the same way as in the usual Bloch case jBBROlj , the Liouville structure al- 
lows to state that the density matrix remains Hermitian through time evolution. 
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Its trace is conserved and it is a positive operator. Hence for a given electric 
field E(t), the elements of the density matrix are bounded, more precisely pop- 
ulations are bounded 

\pUt)\. \p]j{t)\ < Tr(p(t)) = Tr(p(0)), for all i G r and j G I\ 

as well as coherences 

\pm\ < ^jpumAi) ^ all ^, J e /^ 

and likewise 

IpJfcWI, Ip^WI ^ for all » e r and j, fc G r. 

4.2 Coupling with the Maxwell equations 

The density matrix governed by Bloch equations in the previous section was only 
depending on time. We can now consider a collection of quantum dots which 
are scattered in space and interacting, not directly but through the interaction 
with an electromagnetic wave that propagates through the medium. This can 
be modeled by a density matrix, that now depends on time and space, which 
is coupled with Maxwell equations for the laser field through the expression of 
polarization. We therefore address the system 

' ^i^t^^ = — curl e, 

edtE = curlH- 9tP, 

P =A/'bTr(Mp), (17) 

where all the variables depend on time and space in 3 dimensions: the electric 
and magnetic fields E and H, the polarization P, and the density matrix p. Such 
models have already been written and studied mathematically and numerically 
in a few physical contexts. Here the specificity is the fact that the Bloch equation 
is nonlinear in p. Note that even in the case when V does not depend on p, the 
full coupled model is already nonlinear since V is affine in E. 

This system can be cast in the abstract setting of |DS10| . In this paper, 
they introduce a general abstract setting able to treat both Maxwell-Landau- 
Lifschitz and classical Maxwell-Bloch equations. In this setting the electromag- 
netic field is supposed to exist in all space R^. Matter described by the density 
matrix is only occupying a subdomain of M^. The variables are gathered in 
one variable U — (u, v), where u — (tti, M2) = (H, E) and v = p. The variable u 
can be view as 6 real variables and variable v as d = (card(/'^) -I- card(/^))^ real 
variables. This variable U is supposed to be in = L^^K^.j^e^ ^ L'^{n;R'^). 
The abstract system reads 



{dt + B)u ^ • l)F{v, u), for x G M^ 
dtv — F{v, u), for X E i}. 



(18) 
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In this formulation k{x) — {ki{x), K2{x)) = {^i,s). It is uniformly positive as 
needed in |DS10| . The linear differential operator B is defined on i?curi x -ffcuri 
by B{ui,U2) = {ki^ curlM2, — curlui). The variable v is extended by v on 
the whole and is zero outside $1. We can identify li = 0, h = — A/b Tr(M-), 
and F{v,u) = —j^[V{p), p\. The system (17) verifies the hypotheses given in 
[DSTO] . namely 



F is afHnc in u: F{v,u) — Fo{v) + Fi{v)u, 
forj=0,l, F,(0) = 0. 



• for all i? > there exists Cp{R) such that for all v £ Bn (ball of radius 

R in M'*), \Fj{v)\ + \d^Fj{v)\ < Cf{R). 

• There exists K >Q such that for all (m, u) e ^ K'', F{v, u) ■ v < K\u\'^ . 

We have in particular used the L°° bounds of the previous section, more precisely 
we look for 

V £ L°°((Q,oo)-L°"{n-«^)). (19) 
Besides we suppose to have at time t = Q the conditions 

div(KjMj - Ijv) = 0, for j 1, 2, (20) 

which are indeed the physical conditions divB = and div(eE + P) = 
|Dum05) . The structure of Equation (181 ensures that this condition holds 
for all time if it is valid at the initial time. 



4.3 Cauchy problem 

In this section, we state without proof the results obtained in jPSlO] and that 
can apply to our context. The first result addresses the existence of global finite 
energy solutions. 

Theorem 1 (Theorem 3, |DS10j ). For any initial data Uq — (wo,wo) satisfying 
(20), there exists U € C([0,oo);L^) which is a global finite solution to (18)-(20). 



Moreover, for all T > 0, there exists a constant C that only depends on T , F , I 
and ||wo||l°= , such that for all t £ [0, T], ||C/(i)||L2 < C'||C^o||l2 • 

To have a uniqueness result we need some regularity on e and ^. Usually in 
physical contexts e and ^ experiment jumps on the boundary of Vl, but we do 
not know how to tackle with this problem. We hence assume that 

Q. is bounded and - 1 e Cjf (M^), for i = 1, 2, (21) 

which means that is 1 outside 17 and the transition across the boundary is 
smooth. We also assume that the initial data for the electromagnetic wave is 
smooth enough with Sobolev regularity. 

Theorem 2 (Theorem 5, jPSlOp . Under the assumptions of Theorem^ and 



(21), let s e (0, 1] and curluo* e H^-^iM?) for i = 1,2. Then 
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(i) any solution given by Theorem^ satisfies cmlui G C{[0,oo); H" 
for I = 1,2; 



(ii) if s = \, there exists only one solution to (18|-(20l with initial data Uq 



Item (ii) stems from the fact that h = and F does not depend on ui but 
only U2- 

5 Numerical experiments 
5.1 Self-Induced Transparency 

Self-Induced Transparency (SIT) is a typical two-level phenomenon: using a 
light pulse which is resonant with the transition, absorption and stimulated 
emission are combined to obtain exact population inversion and an unchanged 
electric field. This phenomenon has been predicted theoretically and confirmed 
experimentally AE87, MH67, GS70 . 

The propagating field is a pulse given by 

E{t, z) = £{t, z) sin(ajo(i ~ z/v)), 

where v is the velocity of the pulse, ujq is both the center frequency of the pulse 
and the transition frequency of the medium, and £{{t, z) is the pulse envelope. 
It is shown that the envelope is not reshaped by the medium, only if it is a 
symmetric hyperbolic secant 

8{t, z) = £o sech \ — ) , where = ——■ 

\ T J WLT 

In this expression r is the pulse duration and m = M /h, where M is the dipolar 
moment associated to the transition. According to the Area Theorem, the 
medium undergoes k exact inversions if 



/oo 
£{t, z)dt = kiT. 
-oo 



The corresponding pulse is called a fcTr-pulse. For our test-cases, we will use a 
27r-pulse, for which the medium is inversed and goes eventually to its original 
state. The return to the initial state is a easy-to-check criterion to validate 
numerical approaches, as has been already done in |BF06[ [7AG95| . It is easy 
to compute that 

A = mrEf) [arctan(sinht)]'^^ = mrEQir, 
and hence a fcTr-pulse is obtained for the amplitude £q = kh/Mr. 
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5.2 Adaption to the quantum dot context 



In this paper we want in particular to investigate the vahdity of the vanishing 
intra-band assumption. To this aim we need a minimum of three levels and we 
therefore adapt the SIT experiment to our framework. We absolutely do not 
claim that SIT has any practical application in the quantum dot context, but 
only choose this test-case because of the easiness to interpret the results. 




Figure 1: Adaption of the SIT test-case to the quantum dot context, (a) Original 
two-level case; (b) 2 three-level test cases. 

In Figure [ija) the original two-level test case is represented, for which there 
is a single conduction level and a single valence level, separated by the energy 
corresponding to the field frequency. The upper plot represents the (normalized) 
time-evolution of the electric field. The time-evolution of the population of the 
initially empty conduction level is given by the lower curve. We observe that 
the medium undergoes two complete population inversions. 

In Figure [T|b) we have three-level test cases with a single conduction level 
and two valence levels. Coulomb interaction is not taken into account. In the 
first case (represented by solid lines both on the plot and on the scheme) the 
transition between the two valence levels is also resonant with the field and this 
destroys the SIT phenomenon. It suffices to get both valence levels far apart 
enough (e.g. 2wo as in the second case, represented by dashed lines) to recover 
SIT. We use this last configuration as basis test-case for the following numerical 
experiments. 
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5.3 Numerical features 



The simulations are performed using a code based on a finite difference Yee 
scheme and a relevant choice for the time discretization of the Bloch equation 
(see |Bid03| . It allows to keep the good properties of the original Yee scheme: 
second order and explicitness. A splitting scheme, first described in [BBROlj . 
allows to preserve positiveness at the discrete level. It is strongly based on 
the exact solution given by Equation ([s]). It has been adapted to include also 
Coulomb interaction, still preserving positiveness, but at the cost of a loss of 
approximation order, which becomes one. 

Integrating the zero intra-band coherences assumption is a priori a problem 
since it destroys the Liouville structure and an exact solution is no more avail- 
able. It is however possible to solve a Liouville-like equation and set artificially 
intra-band coherences to zero. This adds a step at each time iteration but allows 
to preserve the general structure of the numerical code. 

To determine the right envelope amplitude for numerics, we use the argument 
of |ZAG95j : in practice the input pulse is cut off on an interval t £ [— IOt, lOr], 
therefore the numerical area is 

An = mrSf) [arctan(sinht)]i°^Q = r7iT£^o(0.9999427r), 

which slightly changes the value of Eq. 



5.4 Results 

5.4.1 Impact of Coulomb terms 

To include Coulomb interaction in the 
the coefficients given by Equation ([t]). 
scope of the present paper. We choose 




-R. = 



SIT test case, we have to give values to 
Their exact computation is not in the 
to take them all equal to a single value 




- PS, 



Time (ps) 
(a) Ro = 10-21 



Time (ps) 



(b) Ro = 3.5 10-21 



Figure 2: Impact of Coulomb terms for two interaction strengths. 

For small values of Rq (see the evolution of p^i described for Ro = lO-^i 
in Figure [2ja), solid plot) SIT is only slightly affected. In this figure, the 
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dashed plot corresponds to the reference case {Rq = 0) and is the same as the 
dashed plot of Figure [T|b). The effect is clearer for stronger values of Rq (e.g. 
i?o — 3.5 10~^^ in Figure ^h)). The total inversion is prevented by Coulomb 
interaction. By the time the electromagnetic wave arrives at the space point 
of interest, Coulomb interact has relaxed part of the population of the higher 
valence level P22 to the lowest valence level p[i (solid plot)). Only the remain- 
ing electrons (about half of them in our example) take part in the phenomenon, 
which explain the dashed plot (and also the slight effect in Figure [2]j a)). Al- 
though inversion is not complete, the return to zero of pli is observed in this 
test case. 



5.4.2 Impact of vanishing intra-band coherences 

We first test the impact of the vanishing intra-band coherence assumption on 
the Coulomb-free model. We always use the same experimental setting (see 
Figure [3](a)) and this assumption amounts to taking = 0- 




Figure 3: Impact of vanishing intra-band coherences on the Coulomb- free model. 

The result is displayed in Figure |3jb). The final equilibrium state for matter 
is slightly changed and pl^, which is given by the solid curve, does not eventually 
return to zero. The last inversion is not a total one. 

The test cases with Coulomb terms is much more interesting. If Rq is low 
(see Figure |4][a)), the result is not much affected by this assumption. 
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Figure 4: Impact of vanishing intra-band coherences on the full model. 

In the case of Rq — 3.5 lO-^i , we see in Figure[4|^b) that the final equilibrium 
state is not physical {pli < 0). This is due to the destruction of the Liouville 
structure. We can easily explain why this is not observed when Rq is low. When 
we are close to the SIT experiment we have a typical two-level phenomenon: 
p^j^ ~ during the whole experiment. For a two-level system, the positiveness 
of the density matrix is equivalent to 

• the positiveness of each diagonal term (populations) , 

• the estimation of coherence by populations, here: |p52^P — Pii/°22' 
second valence level would be the only relevant one (see |BBR 01 ). 

Setting intra-band coherences to zero within the numerical process does not 
affect these properties, and the iteration use in the proof of the positiveness of 
the density matrix applies. 

But for a three-level system (and the case Rq — 3.5 lO-^i is a true three-level 
case) the positiveness of the matrix involves some more properties, which are 
affected by setting to zero. Although trace is still conserved, the positiveness 
of the population is affected. We would of course have the same result with a 
dedicated code where intra-band coherences would simply not be computed. 
Besides the effect is clear enough in Figure |4][b) not to be attributed to simple 
round-off errors. The conclusion is that even if intra-band coherences seem not 
to be very relevant for some physical applications, it is very important to include 
them in the mathematical description and in the numerical computation to keep 
the natural mathematical structure of the density matrix. 

Remark 2. In absence of electromagnetic field the evolution equation is reduced 
to 

ihdtpli = 0. 

This is only due to the fact that there is only one conduction level in our test 
case, and does not depend on the intra-band coherence vanishing assumption, 
nor on the fact that we took a single Rq Coulomb coefficient. Hence, when the 
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population of the coherence level has been set into a non-physical state (and 
this is due to the intra-band coherence vanishing assumption), it remains in this 
state for ever. 

6 Conclusion 

In this paper, Bloch-type equations have been derived considering Coulomb 

effects in quantum dots. We have shown analytically and numerically that 
Coulomb effects are not negligible in some quantum dot structures, and we have 
given the link between mathematical properties and physical relevancy of the 
Bloch model and more specifically in the treatment of intra-band coherences. 
Then this model has been coupled with the description of laser propagation 
in the quantum dot structures, leading to a Maxwell-Bloch system for which 
we have studied the local and global Cauchy problems. This system has been 
implemented numerically and simulations have been performed on a self-induced 
transparency test-case. In particular, we have tested the impact of Coulomb 
parameters and intra-band coherences. Further work will include additional 
effects in the same Bloch-type framework. 
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A Canonical form for ^ 



ci\ *J^2 «! Q-^ aia^ 



-Sa2a['"aivl, + Was ^ai 



+'^ala^U^'^^'a2 - ^^I; ?^ai 



B Commutator of cic, with iif'^ ^ 



C Commutator of v^^Ci with i^'^ ^ 



Sja'^ 4i Ci + 4^ Vc,'^ v] C^^ ( 
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